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Due to strong interlayer correlations, the bilayer quantum Hall system is a single coherent system 
as a whole rather than a weakly-coupled set of two independent systems, which makes conventional 
tunneling theories inapplicable. In this paper, we develop a theory of interlayer tunneling in coherent 
exciton condensates of bilayer quantum Hall systems at total filling factor vt = 1. One of the most 
important consequences of our theory is that the zero-bias interlayer tunneling conductance peak 
is strongly enhanced, but fundamentally finite even at zero temperature. We explicitly compute 
the height of the conductance peak as a function of interlayer distance, which is compared with 
experiment. It is emphasized that the interlayer distance dependence of the conductance peak 
is one of the key properties distinguishing between the spontaneous coherence due to many-body 
effects of the Coulomb interaction and the induced coherence due to the single-particle tunneling 
gap. It is also emphasized that, though the strongly enhanced tunneling conductance originates from 
the interlayer phase coherence, it is not the usual Josephson effect. We propose an experimental 
setup for the true Josephson effect in couterflowing current measurements for a coupled set of two 
bilayer quantum Hall systems, which is a more precise analogy with the real Josephson effect in 
superconductivity. 

PACS numbers: 73.43.-f, 73.21.-b 



I. INTRODUCTION 

Strongly-enhanced interlayer conductance peaks ob- 
served by Spielman et alX near zero bias in bilayer quan- 
tum Hall systems at total filling factor Vt = 1 have at- 
tracted much interest for many reasons. One of the main 
reasons is the existence of a rather precise analogy be- 
tween the bilayer quantum Hall effect 2 and superconduc- 
tivity. That is, the ground state of the bilayer quan- 
tum Hall effect at small interlayer distance d/ls <C 1 
(Ib = y/hc/eB) can be mapped onto the BCS-type wave- 
function describing the Bose-Einstein condensate of ex- 
citon pairs formed between electrons and holes residing 
across the interlayer barrier. 

Bosc-Einstcin condensation of neutral excitons has, in 
fact, been sought after in semiconductors for decades. 
In particular, there have been fascinating recent exper- 
iments on possible condensation of optically generated 
indirect excitons^, for which, however, evidence is not 
yet conclusive. On the other hand, in the bilayer quan- 
tum Hall effect, it is generally accepted that the strongly- 
enhanced conductance peak is a direct indication of the 
macroscopic phase coherence. 

To be concrete with respect to the mapping between 
superconductivity and the bilayer quantum Hall effect, 
let us begin by writing the ground state wave function 
at d/ls — * 0. The explicit ground state wave function 
at (L/Ib —> is very instructive in illustrating the essen- 
tial physics even though the ground state is substantially 
more complicated for general d/ls, especially near the 
critical d/ls where the exciton condensate disappears. 
(Note that we will use numerical methods to study the 



ground state at general d/ls-) The ground state wave 
function at d/ls — > is known as the Halperin's (1,1,1) 
stated: 

i^ii>=n( c ™T+ c ™i)i°>' 

ra 

where m is a momentum index of the lowest Landau level. 
In the above, the pseudospin representation is used: f 
and I indicate the top and the bottom layer, respectively. 
It is important to note that, contrary to the usual repre- 
sentation of the (1,1,1) state, 

ipiu = II ( z * ^ z i) \\ ( z k - zi) Yl ( Zm ~ z ") ( 2 ) 

i.jeT Me| me\,nel 

which denotes only orbital components, Eq.JD describes 
the full wave function including both the orbital and the 
layer degree of freedom^. 

The ground state wave function described by Eq.lfT|) 
has an isomorphic structure to the BCS wave function, 
which can be made more apparent after the following 
reorganization: 

i^m) = n( i+c -T c «u)rK'j°> 

m m' 

= JI(1 + c I t nT c mi )|new vacuum) 

m 

= 11^ + c mT /l Ll)l new vacuum), (3) 

m 

where h) is a creation operator for holes, acting on the 
fully-filled bottom layer. Because of this analogy, it is 



rather natural to expect that the bilayer quantum Hall 
state may have a coherence in the phase associated with 
interlayer electron-number difference, which is similar to 
the phase coherence associated with the total number 
of Cooper pairs in superconductivity. The phase coher- 
ence between states with different Cooper-pair numbers 
is the origin of the Josephson effect in superconductivity. 
Naturally, this similarity led many previous authors^ 
to predict the Josephson effect in bilayer quantum Hall 
systems. So, in this context, the strongly enhanced con- 
ductance observed by Spielman et alX seemed to be ex- 
actly the experimental verification needed. There are, 
however, some important properties of the conductance 
peak indicating that this phenomenon is not the con- 
ventional Josephson effect; most notably, experiments 
suggest that the height and the width of the zero-bias 
conductance peaks saturate to finite values when avail- 
able finite- temperature data are extrapolated to the zero 
temperature limilj^. This means that there is no DC 
current strictly at zero bias voltage in contrast to the 
real Josephson effect in superconductivity. 

This apparent discrepancy gave rise to two groups of 
thought. In one group, the enhanced conductance is still 
regarded as DC Josephson effect, but its height is reduced 
by complicated disorder-induced fluctuationsi2*ii*i2ii^. 
On the other hand, others^ have argued that there is no 
exact analog of the Josephson effect in the experimental 
setup measuring interlayer tunneling currents in a sin- 
gle set of bilayer quantum Hall systems. It is because 
the bilayer system as a whole is a single Bose-Einstein 
condensate(BEC), not a set of two independent BEC's. 
While this argument itself is clearly true, it is still nec- 
essary to explain why and how the enhanced interlayer 
conductance is related to the BEC of excitons. The res- 
olution to this crucial issue has remained elusive even 
after a very extensive body of research efforts continuing 
to recent yeara 15 i 16 i 17 i 18 i 19 i 20 . 

Important issues to be addressed are, in particular, (i) 
exactly what is the process of coherent tunneling in inter- 
layer tunneling current measurements, (ii) how this pro- 
cess can be mathematically formulated, and finally (iii) 
what is the precise relationship between the interlayer 
phase coherence and the zero-bias conductance. It is the 
goal of this paper— to address these questions. Specif- 
ically, we would like to provide a quantitative analysis 
of the zero-bias conductance peak as a function of cI/Ib, 
which, in turn, generates a sharp distinction between the 
coherent and incoherent tunneling processes. 

The d/ls dependence of the conductance peak is im- 
portant also because it distinguishes between the coher- 
ence due to many-body effects of the Coulomb interac- 
tion and that due to the single-particle tunneling. In this 
paper, we are primarily interested in the regime of suffi- 
ciently small single-particle tunneling, i/(e 2 /e?s) -C 1. 
We are interested in this regime mainly because ex- 
perimental values of tj '(e 2 / 'tis) are really quite small; 
in recent experimental setups, it can be as small as 
10~ 6 — 10~ 7 . This limit is also quite interesting in a theo- 
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FIG. 1: Diagram showing the difference between (a) coherent 
tunneling and (b) incoherent tunneling. Note that, for the 
coherent ground state, there is no chemical potential differ- 
ence between the two layers, as depicted in (a). When d is 
sufficiently large as shown in (b), however, each layer becomes 
independent, and therefore the interlayer coherence is lost. 

retical point of view since the interlayer phase coherence 
is induced purely by many-body effects of the Coulomb 
interaction in the absence of single-particle tunneling. In 
later sections, we will investigate similarities and differ- 
ences between the spontaneous coherence caused by the 
Coulomb interaction and the induced coherence caused 
by single-particle tunneling. 

This paper is organized as follows; we begin in Sec- 
tion [H] by distinguishing physical situations for coherent 
tunneling from those for incoherent one. We provide a 
physical picture as to exactly how interlayer tunneling 
can be enhanced in the presence of coherence. We, then, 
present our quantitative analysis in Section II I II by pro- 
viding a mathematical derivation of the many-body tun- 
neling Hamiltonian and its connection to the interlayer 
tunneling conductance. In Section IIVI we numerically 
evaluate the coherent tunneling conductance as a func- 
tion of d/ls and compare it with experiment. In Section 
fVl we propose an experimental setup for the Josepson 
effect in counterflowing current measurements for a cou- 
pled set of two bilayer quantum Hall systems. We finally 
conclude in Sec. I VII 

II. COHERENT TUNNELING VERSUS 
INCOHERENT TUNNELING 

As mentioned in the introduction, the bilayer quan- 
tum Hall system itself is a single superfluid system. The 
very existence of the excitation gap required for the quan- 
tized Hall resistance indicates that the ground state wave 
function is robust against small perturbations, in which 
case the bilayer system can be assumed to maintain its 
equilibrium. External perturbations such as bias voltage, 
then, can be taken as perturbations to the bilayer system 



as a whole, not to the individual layers. In this situation, 
there is no chemical potential difference between the two 
layers, as depicted in FigE](a). 

In the situation described by Fig|T] (a), tunneling is 
coherent because the aforementioned robustness of the 
ground state is due to the interlayer coherence. Physi- 
cally speaking, coherent interlayer tunneling is obtained 
as follows: Imagine that an electron is inserted into the 
top layer of the bilayer quantum Hall system, as depicted 
in the top diagram in FigEl The inserted electron, then, 
becomes a quasiparticle on top of the ground state with 
an energy penalty (the second diagram from the top in 
Fig HI- The energy penalty can be caused either by the 
Coulomb interaction (for the spontaneous coherence) or 
by the single-particle tunneling gap (for the induced co- 
herence). The so-introduced quasiparticle, then, immedi- 
ately becomes a part of the coherent ground state by res- 
onating between the two layers (the third diagram from 
the top). The electron is finally ejected from the bot- 
tom layer, and thereby the interlayer current flows. It is 
important to note that, when there is an interlayer co- 
herence, tunneling can be strongly enhanced because all 
inserted electrons will tunnel in a synchronized fashion. 

The situation is quite different when there is no inter- 
layer coherence. There are two important cases where 
the interlayer coherence can be lost; (i) for a fixed filling 
factor, namely vt — 1 in our case, cL/Ib becomes large or 
(ii) simply B = in which interactions play a weaker role. 
For the first case, each layer forms an individual compos- 
ite Fermi sea 22 i 23 i 24 at filling factor v — 1/2. What actu- 
ally tunnel between the layers, however, are real electrons 
rather than composite fermions. In this situation, the in- 
terlayer tunneling is suppressed due to the fact that the 
composite Fermi sea is a highly correlated state so that 
the sudden insertion of an uncorrelated electron requires 
a large energy cost. 

When B = 0, the situation is simpler: ordinary elec- 
trons tunnel. In this situation, the interlayer tunneling 
current can be evaluated in a conventional perturbation 
theory assuming that the single-particle tunneling gap, 
t, is small and therefore the two layers become indepen- 
dent and have individual chemical potentials, as depicted 
in Fig|2](b). To be concrete, let us start with the tun- 
neling Hamiltonian, H t : 



(4) 



where t is the single-particle tunneling gap and the pseu- 
dospin representation is used. It is assumed that mo- 
menta parallel to the two-dimensional plane are con- 
served; in other words, there is no tunneling between 
different k's. In this representation, it is not too difficult 
to show that the tunneling current operator, J, is given 
as follows: 
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When t is small, one can use a conventional first-order 
S-matrix expansion to compute the expectation value of 
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FIG. 2: Schematic diagrams for coherent interlayer tunneling. 
Note that each ellipse indicates the resonance between the top 
and the bottom layer. Imagine that an electron is inserted 
into the top layer of the bilayer quantum Hall system, as 
depicted in the top diagram. The inserted electron, then, 
becomes a quasiparticle on top of the ground state with an 
energy penalty (the second diagram from the top). The so- 
introduced quasiparticle, then, immediately becomes a part 
of the coherent ground state by resonating between the two 
layers (the third diagram from the top). It is important to 
note that the interlayer resonance can be caused either by 
many-body effects of the Coulomb interaction or by the single- 
particle tunneling. Finally, the electron is ejected from the 
bottom layer, and thereby the interlayer current flows (the 
bottom diagram). 



the current operator: 

1(f) = -i f dt'([J(t),H T (t')}}, (6) 

which, after some algebra, leads to the following expres- 
sion: 

I = 2et 2 / ^ E m(eKi(e + eV)[f(s) - f(e + eV)]{7) 
J k 

where pw(s) is the density of states for a given k (which is 
also known the spectral function), and /(e) is the usual 
Fermi-Dirac distribution function. Note that PkT( £ ) — 
/9kj.(e) = Pk(£) because of symmetry. Now, since the z- 
component motion of electrons is frozen due to the con- 
finement to the two-dimensional plane, /Ck(e) is sharply 
peaked around its center (which depends on k). As a 
consequence, the integrand in Eq.J7|) vanishes except for 
small eV: 



Pk(e)pk(e + eV r ) ~ pl(s)e 



-(eV) 2 /2S 2 



(8) 
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where the size of 5 is determined by the static disorder, 
which controls the width of the quasiparticlc coherence 
peak. For small <5, as is expected in high-mobility 2D 
systems of interest, this single-particle tunneling peak 
is almost a delta-function with a small resonance width 
eV ~ S. 

This resonance behavior of the interlayer tunneling 
conductance is, in fact, observed at B — in the same 
sample 1 showing the enhancement of the zero-bias tun- 
neling conductance at vt — 1 for small d/ls- So, at least 
superficially, incoherent tunneling seems to generate sim- 
ilar physical consequences as coherent tunneling. There 
are, however, indications suggesting that the B = con- 
ductance peak is actually different from the = 1 peak. 
One indication is that the B = conductance peak is 100 
times smaller than the peak at vt — 1. Furthermore, the 
-6 = peak is at zero bias only when the two layers 
have the same density. If the densities are different, the 
B = peak shifts accordingly. In contrast, the Vt = 1 
peak is locked on to V — 0, even for (small) non-zero 
density differences. Despite these indications, however, 
it is very important to know whether coherent tunneling 
can give rise to results which are manifestly physically 
distinct from the corresponding incoherent tunneling re- 
sults, and, if so, what those results are. 

One of such distinctions can be obtained in the d/ls 
dependence of the zero-bias tunneling conductance. In 
Eq.0, the only dependence of the incoherent tunneling 
current on d is through the tunneling gap, t. Experi- 
mentally, d/ls is varied usually by changing the mag- 
netic field for the same physical sample (Note that the 
electron density is also adjusted to maintain the fixed 
filling factor). In this experimental setup, t should not 
change much since t depends on the physical distance d, 
not on d/ls- Experimentally observed, however, is an 
abrupt change in the zero-bias tunneling conductance: 
the zero-bias conductance peak completely disappears af- 
ter a certain critical distance d c jl^. This is the reason 
why it is so important to develop a quantitative theory 
to compute the zero-bias conductance peak as a function 
of d/ls- We provide such a theory in the next section. 
We believe that the key to understanding the interlayer 
coherent tunneling experiment of Spiclman et aZi is its 
d/lB dependence, which, as a matter of principle, can 
distinguish between incoherent single-particle tunneling 
and the many-body correlation induced coherent effect 
of bilayer excitonic condensates. 



III. DERIVATION OF THE MANY-BODY 
TUNNELING HAMILTONIAN 

In this section, we develop a theory of coherent tun- 
neling which occurs when there is an interlayer coher- 
ence in the ground state. An important starting point is 
that the two layers cannot be regarded as independent 
systems. Instead, they must be treated as a coherent 
whole. Since there is no interlayer chemical potential in 



a single coherent bilayer system, there is no electromotive 
force within the bilayer system and therefore any current 
should be induced from outside. So, it is necessary to 
take into account external leads, as schematically shown 
in Fig^ This, of course, makes any quantitative pre- 
diction dependent on the way in which bilayer systems 
are connected to external leads. However, it is still pos- 
sible to make quantitative predictions on some essential 
aspects of coherent interlayer tunneling. In particular, 
we study the d/ls dependence of the zero-bias tunneling 
conductance peak. In doing so, we also show that the 
width of the conductance peak is fundamentally finite 
even at zero temperature, and it is controlled ultimately 
by very small, but finite single-particle interlayer tunnel- 
ing gap t. This makes sense since a finite t is essential to 
experimentally drive an interlayer current. 

Let us begin our quantitative analysis by writing the 
total Hamiltonian including (i) the interlayer (single- 
particle) tunneling Hamiltonian, H t , (ii) the Hamiltonian 
for the Coulomb interaction between electrons, -Hcoui, 
(iii) the Hamiltonian describing the left and the right 
lead, Hl and Hr respectively, and finally (iv) the Hamil- 
tonian for tunneling between leads and the bilayer sys- 
tem, H': 

H = H + H' + H R + H L , (9) 
H = H t +H Couh (10) 

H t = t^ic^Cjnl + c'^Crn^) = 2tS x , (11) 
m 

H' = ]rr KT (fc,m)[ c y£0 CmT +H.c.] 

k , m 

+ ^T Li (p,m')[c[(p)c ma +H.c.], (13) 

p,m' 

where, as before, the pseudospin representation is used 
and m is a momentum index of the lowest Landau level. 
Vhhh is the lowest Landau level projection operator. 
Tji\{k,m) is the tunneling amplitude between the state 
with momentum k in the right lead, and the state with 
to in the top layer of the bilayer system. T^^ip, m') is 
similarly defined. Hr and Hl describe external leads 
as normal Fermi liquids. Note that H t is related to the 
pseudospin magneization in the x-direction, S x : 

Sx= 2 ^2( C U Cm i + c Li c "T)- ( 14 ) 

Also, note that the internal interlayer current operator is 
related to the pseudospin magneization in the y-direction, 
S y : 

./inter = e ti^2( c i^ c ml ~ C^Cmj) EE 2etS y . (15) 

m 
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Now that we are interested in the limit of zero 
t/(e 2 /els), we investigate the possibility of a current op- 
erator which survives as t/(e 2 /els) — * 0. This operator 
will, then, measure the spontaneous interlayer current 
as opposed to the single-particle hopping induced cur- 
rent. To this end, let us consider the following tunneling 
Hamiltonian: 



T Eg — H — Hr — Hl 



(16) 



where E g is the ground state energy of H + Hr + Hl . 
The idea is that, by adding an electron to the top layer 
and removing another from the bottom layer, H' T de- 
scribes the total current flowing through the bilayer sys- 
tem. Note that H' T is second order in H 1 since it is the 
lowest order that can carry the current from one lead to 
the other lead in the limit of zero t/(e 2 I (1b)- 

Normally, the above process would not cause actual in- 
terlayer tunneling, but instead it would result in charge 
build-up. Our bilayer quantum Hall system, however, 
is special such that it has spontaneous interlayer coher- 
ence for small interlayer distances and therefore electrons 
can move back and forth between the layers even in the 
limit of zero t/(e 2 /els)- The total interlayer current is, of 
course, zero in the ground state since the back and forth 
currents cancel each other. Non-zero interlayer currents, 
however, can be induced by applying a bias voltage which 
breaks the balance between the back and forth motion in 
equilibrium. 

Once interlayer currents begin to flow as described by 
the process oiH' T in Ea. 1(16(1 . the currents can flow steady 
only when it is carried by the internal interlayer current: 
I = (>/intcr), as dictated by the charge conservation. The 
exciton condensate can carry the interlayer current by 
adjusting its interlayer phase difference ef>. To under- 
stand this, it is instructive to consider the ground state 
wavefunction at d/ls — * in the presence of non-zero 
interlayer phase: 



l^m(0)>=n( c U+ e ^i)l°>- 



(17) 



It can, then, be shown that the ground state in Ea. f!7p 
carries a non-zero interlayer current: 



(Winter) = 2et(S y ) = 2et(S) sin0, 



(18) 



where (S) = \{T, m 4nt c ml)\- _ 

It is important that there is a maximum current that 
can be carried by the coherent state: 



(Winter) "^critical — '2(',t\S). 



(19) 



The existence of the maximum coherent current is analo- 
gous to that of the critical current in the superconductiv- 
ity. This termination of the coherent tunneling current 
beyond a critical value is the reason why the tunneling 
conductance has a narrow peak near zero bias. That 
is, the tunneling current is coherent only within a small 



window of bias voltage. Once the bias voltage gets larger 
than a critical value, the tunneling current becomes too 
large to be carried only by coherent tunneling. The ex- 
cess current must, then, be carried by incoherent tunnel- 
ing, which has a low conductance, as argued previously 
in terms of the formation of the composite fermion sea. 

It is important to note that H' T in Eq.(16) is zeroth 
order in t whereas, in usual incoherent tunneling, the ef- 
fect of the tunneling Hamiltonian would be second order 
in t since the expectation value of S x in Ea. 1)11(1 is also 
proportional to t. Therefore, for spontaneously coherent 
systems, the interlayer tunneling conductance does not 
depend on t in the limit of zero t while the charge con- 
servation condition simply imposes an upper bound on 
the value of coherent tunneling current, which is propor- 
tional to t. In other words, our theory is consistently 
constructed up to the first order of t, which is contrasted 
to the second order behavior of the usual theory for in- 
coherent tunneling conductance. 

Now, to proceed further, let us examine H' T in Ea. 1(16(1 
in detail. Since the bilayer quantum Hall state is in- 
compressible at sufficiently small d/ls, adding or remov- 
ing electrons costs a finite energy, Aj2S which is equal 
to either (i) the Coulomb self-energy of a quasiparti- 
cle, Ac, for the spontaneous coherence occurring when 
t <C e 2 /elB, or (ii) the single-particle tunneling gap, t, for 
the induced coherence occurring when t » e 2 /els- While 
it is possible to investigate both regimes, the spontaneous 
coherence is particularly interesting, as seen in later sec- 
tions. It is, however, sufficient at this stage to know that, 
no matter whether A is Ac or t, A is independent of the 
momentum index, m, in the lowest Landau level. So we 
just replace H + Hr + Hl — E g by A. 

Note that we are able to make the above simplification 
because H' in Ea. 1(13(1 is chosen such that electrons are 
inserted directly into Landau level eigenstates with mo- 
mentum indices m. Since each Landau-level eigenstate 
remains to be a well-defined quasiparticle eigenstate in 
the presence of the Coulomb interaction, one may take 
each m state to be a very good approximation to the ex- 
act eigenstates of H + Hr + Hl with eigenenergy E g + A 
(Note that A is a self-energy correction). 

Now, we assume that the tunneling amplitudes 
Tfi-f (k, m) and Tli(j>, m') are more or less independent of 
momenta k and p, which is a common practice in tunnel- 
ing theories where tunneling occurs only within a narrow 
region of energy near Fermi surface. Keeping only the 
terms in H' T relevant for transporting electrons from one 
lead to the other, we arrive at the following many-body 
tunneling Hamiltonian: 



Ht — 



J2[cU k )cL( P )f^ +cl(p)c R (k)f 



where 



t =a^ t - 



(20) 



(21) 



d/l 



FIG. 3: Feynman diagram for coherent tunneling in bilayer 
quantum Hall systems. The vertex operator T contains many- 
body effects of the exciton condensate. Trl is the tunneling 
amplitude and Ac is the Coulomb self-energy of a quasipar- 
ticle. 



and Trl(tti) = TR^(kp, m)Tn(kF, to). Based on Ht, 
the tunneling current operator J is defined as follows: 



J = ei 



[4(fe)ci(p)Tt _ cl( p ) CR ( k )f] . (22) 



Assuming that the tunneling amplitude between leads 
and the bilayer system is small, one can compute the 
expectation value of current operator via a conventional 
first-order S-matrix expansion: 



I(t) = -i dt'([J(t),H T (t')}). 



(23) 



The new aspect of our tunneling theory is the vertex 
operator, T, which contains all of many-body effects of 
the exciton condensate. Ea. (|23[l can be evaluated further 
using the Feynman diagram depicted in Fig|3 

/°° dp 
^-A R (k,e)A L (p,e + eV) 



x lf(e)-f(e + eV)} 
= 4ne 2 D R D L \{f)\ 2 V 



(24) 



where Ar (A^) is the spectral function of the right (left) 
lead, /(e) is the usual Fermi-Dirac distribution function, 
and Dr (Dl) is the density of states at the Fermi sur- 
face of right (left) lead. Note that we have used \(T)\ 2 
in Eq.lPU instead of {T^T) since |(T)| 2 is a much more 
acute measure of the spontaneous interlayer order in 
finite-size system diagonalization studies while the two 
quantities generate identical results in the thermody- 
namic limit. 

Ea. l24|) indicates that there is no DC Josephson ef- 
fect in a conventional sense since the conductance G 
(= dl/dV oc |(T}| 2 ) is finite. It is important to know, 
however, that the coherent tunneling current is zero if 
(T) = 0. Remember that (T) measures the phase co- 
herence between states with various interlayer number 



FIG. 4: Interlayer coherence quantified by the pseudospin 
magnetization, (S x ), for the AT = 11 system. (S x ) is com- 
puted via exact diagonalization as a function of the interlayer 
distance, djls, and the single-particle tunneling parameter, 
t/(e 2 /eIb)- Note that the shape of (S x ) is qualitatively similar 
to the usual magnetization curve as a function of temperature 
and external magnetic field, which correspond to interlayer 
distance, d, and single-particle tunneling gap, t, respectively. 
Also, note that {S x } in the plot shows typical finite-system 
smearing of a sharp phase transition. The thermodynamic 
limit is estimated in Fig|S|for t/(e 2 /cIb) = 0. 



differences, N Ie i (Note that T oc c^Cmi). So, unless the 
ground state is a coherent linear combination of states 
with various Ar re i, (T) is zero, and so is the tunneling 
current. As mentioned before, this is precisely analogous 
to the phase coherence between different number eigen- 
states in superconductivity, which is responsible for the 
Josephson effect. In this sense, the interlayer tunneling 
conductance is related to the Josephson effect. However, 
we emphasize that the conductance is fundamentally fi- 
nite even at zero temperature. In the next section, we will 
actually evaluate the interlayer tunneling conductance as 
a function of d/ls by using numerical methods. In partic- 
ular, we will be interested in the normalized conductance 
since the absolute scale of the conductance is sensitive to 
sample-specific details such as Dr, Dl and Trl- 

In concluding this section, we mention that the ba- 
sic idea in our derivation of the many-body tunneling 
Hamiltonian is rather similar to that of the linear re- 
sponse theory; if a perturbation is small enough, one 
can extract various dynamic properties of the biased sys- 
tem from its ground state properties in equilibrium. In 
our manuscript, we make a connection between the pseu- 
dospin order parameter of the ground state and the tun- 
neling conductance of the steady state. This connection 
can be valid as long as the interlayer current is not too 
large. 



IV. NUMERICAL EVALUATION OF THE 
COHERENT TUNNELING CURRENT 



To study the d/ls dependence of the tunneling con- 
ductance, one has to compute |(T)| 2 in Ea.l|24|l. which 
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can be further reduced as follows 



A 2 



(25) 



where N is the total number of electrons. Here, we 
have simplified the evaluation of the many-body tunnel- 
ing matix element, T, by assuming that the interlayer 
phase coherence is uniformly obtained throughout the 
system; in other words, the interlayer coherence order 
parameter {c\ n ^c m {) is independent of m. 

In this situation, the m-dependence of tunneling am- 
plitudes between leads and correspoding adjacent layers, 
Trl(ti), affects only the overall magnitude of the tunnel- 
ing conductance, leaving its d/lg dependence unchanged. 
Note that effectively what we are assuming is that the 
lead does not alter important bilayer properties such as 
the interlayer distance and the single-particle tunneling 
amplitude. Now, since (c^Cmi) is independent of m, 



) = (S x )/N. S x [= E m (cLc m i + c 



lJ. C ™T 



)/2] 



can be regarded as an order parameter of the exciton 
condensate. It is natural to assume that ^ m T RL (m) /N 
does not depend on cI/Ib since Trl is only dependent 
on the hopping amplitude between lowest-Landau-level 
states in a given layer and plane-wave states in the ex- 
ternal lead connected to the layer. The d/ls dependence 
of the conductance is, therefore, solely determined by 
(5,) 2 /A 2 . 

Figure 0] shows the three-dimensional plot of (S x ) as a 
function of the interlayer distance, d/ls, and the single- 
particle tunneling gap, t/ (e 2 / cIb)- {S x ) in the plot was 
computed via exact diagonalization for the N = 11 sys- 
tem. Note that, when computing (S x ) in finite systems, 
it is very important to take into account fundamental 
fluctuations in N le \; the true ground state is a coherent, 
linear combination of states with various -/V re i. Technical 
details can be found in the literature 2 ^^* 2 ^. Also, note 
that (S x ) in the plot shows typical finite-system smear- 
ing of a sharp phase transition. It is interesting that 
the shape of {S x ) is qualitatively similar to the typical 
magnetization curve as a function of temperature and 
external magnetic field, which respectively correspond to 
d/ls and t/ (e 2 /els)- 

Figures |H and show that, for d < dc — 1.4 — 1.7Zb, 
{S x ) does not vanish even in the limit of t/(e 2 /els) — > 0, 
which is, by definition, the evidence for the spontaneous 
order. Since (5a;) becomes zero for d > dc in the thermo- 
dynamic limit, it is natural that the tunneling conduc- 
tance becomes zero as well, as dictated by Eas. (|24H and 
(|25p. On the other hand, if induced by the single-particle 
tunneling alone (or, if t ^> e 2 /e^), the coherence does 
not completely vanish at a finite d/lg, nor does the cor- 
responding coherent tunneling conductance peak. This 
is not consistent with experiment where the conductance 
peak disappears at a finite critical interlayer distance. 

So, from now on, we focus on the limit of very small 
single-particle tunneling gap, i/(e 2 /eZs) — > 0, where the 
coherence is spontaneous. Fig[51plots {S x ) 2 as a function 




FIG. 5: Normalized expectation value of the condensate order 
parameter {S x ) 2 in the limit of t/(e 2 eZ_g) — > 0. In order to 
incorporate finite-system corrections, {S x ) is normalized in 
such a way that it is divided by (N + l)/2 for TV odd and 
by \/N(N + 2)/2 for N everA N is the total number of 
electrons in finite-system exact diagonalization studies. The 
shaded region indicates the thermodynamic-limit estimate for 
(&> 2 - 
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FIG. 6: Energy cost of adding (removing) an electron into 
(from) the ground state of the bilayer quantum Hall system. 
Note that this energy cost is essentially the self-energy of a 
quasiparticle due to the Coulomb interaction. The Coulomb 
self-energy is computed via exact diagonalization (for various 
finite systems) as a function of interlayer distance, d/ls- For 
comparison, also plotted is the self-energy estimate obtained 
in the Hartree-Fock approximation. 



of d/ls for various particle numbers when tj '(e 2 / els) 
0. While the accurate estimate of the thermodynamic 
limit of {S x ) 2 is difficult due to numerical uncertainties, 
it is reasonable to assume that the true thermodynamic 
limit lies within the shaded region in Fig[SJ As typical in 
finite-system calculations, inflection points are taken to 
be the finite-system signature for the true critical point 
in the thermodynamic limit. 

In FigH3 we compute the Coulomb energy cost of cre- 
ating a quasiparticle, Ap, which is the only energy cost 
when t/(e 2 /el B ) -> 0. That is, A = A c in Eq.lj^. A c 
is computed via exact diagonalization in finite systems. 
It is interesting to note that, for general values of d/ls, 
the exact Ac is much lower than the estimate obtained 
in the Hartree-Fock approximation 29,30 . 

Finally, in FigEI our theoretical estimate for the nor- 
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FIG. 7: Normalized interlayer tunneling conductance peak 
(the shaded region in the plot) as a function of d/la in com- 
parison with experimental data from Spielman et alp*. We de- 
fine the normalized conductance as the conductance divided 
by its maximum value as a function of d/ls- The shaded re- 
gion is obtained from the thermodynamic estimate of (S x ) 2 
in FigEl 



malized interlayer tunneling conductance near zero bias, 
i.e. (S x ) 2 /Ac, is compared with experimental data of 
Spielman et alX. We define the normalized conductance 
as the conductance divided by its maximum value as a 
function of d/lg- The shaded region in the plot is the 
thermodynamic estimate. To the best of our knowledge, 
the direct comparison with experimental data for the con- 
ductance peak is made here for the first time. It is inter- 
esting that our theory predicts that the conductance peak 
decreases as d/ls decreases below roughly 1.2. Remem- 
ber that the decrease in the conductance peak at small 
d/ls is due to the increase in Ac while the pseudospin 
magnetization is saturated. 

We would like to emphasize that our theoretical esti- 
mate of the conductance peak has only one fitting param- 
eter: the overall prefactor dependent on sample specifics. 
Considering that the only fitting parameter is an overall 
scale factor, we argue that the agreement with experi- 
ment is not too bad. Specifically, our theory provides 
a reasonably accurate estimate for the critical distance 
at which the conductance peak disappears. While the 
critical distance was estimated previously, its determi- 
nation has been based on indirect evidence: it was de- 
termined either (i) by studying the collapse of the over- 
lap between the Halperin's (1,1,1) state and the exact 
ground stat o 31 ' 32 , or (ii) by studying the collapse of the 
low-energy excitation in a time-dependent Hartree-Fock 
approximation 3 ^ or (iii) by approaching the critical point 
from the incoherent side and studying the effect of the 
single-particle tunneling 33 . Our theory, on the other 
hand, is a direct study of the tunneling conductance it- 
self, as measured experimentally, rather than a consid- 
eration of an energy gap collapse or similar purely the- 
oretical constructs associated with the quantum phase 
diagram. 

Aside from the direct estimation of the critical dis- 
tance, our theory predicts a rather steep rise of the peak 



height as djlu is reduced from the critical value, as seen 
in Fig. [7| This steep rise cannot be explained by the 
collapse of the pseudospin magnetization, (S), alone; it 
requires a detailed consideration of the coherent tunnel- 
ing process. 

In the preceding section, we have shown that there is 
a critical interlayer current allowed without breaking the 
phase coherence. Consequently, for a sufficiently large 
bias voltage, the coherent interlayer current should be 
cut off. We have argued that this is the reason why the 
interlayer conductance has a shape of the sharp peak near 
zero bias. Now, we would like to discuss what determines 
the width of the conductance peak. By equating Eas. ijT§|) 
and H'24|) ■ one can show that, while the proportionality 
constant strongly depends on sample-specific details, the 
width of conductance peak is proportional to t. Due to 
the strong dependence on sample details, the accurate 
estimation of the critical current is beyond the scope of 
this paper. It is, however, encouraging to find that typ- 
ical width of the conductance peak (~ lO^eF) is in the 
similar order as the single-particle tunneling gapAi Our 
finding that the width of the conductance peak is propor- 
tional to the tunneling amplitude t shows that the exis- 
tence of a finite symmetry breaking term, i.e. a nonzero t, 
is essential for the experimental observations in Spielman 
et al. 1 



V. PROPOSAL FOR THE JOSEPHSON EFFECT 
IN CONTERFLOWING CURRENT 
MEASUREMENTS 

Until now, we have studied interlayer tunneling in a 
single bilayer system, which, we showed, is not an exact 
analog of the Josephson effect. In this section we propose 
a much more direct analog of the real Josephson effect. 
Since the coherent particle-number fluctuation in super- 
conductivity holds a parallel with the coherent fluctua- 
tion of interlayer number differences in bilayer quantum 
Hall systems, it is natural to expect that the correspond- 
ing Josephson effect for the bilayer quantum Hall systems 
must manifest itself in current measurements associated 
with interlayer number difference. Such measurements 
are counterflowing current measurements 34,35 . 

To be specific, let us consider a weakly coupled set of 
two bilayer quantum Hall systems, say A and B (four lay- 
ers altogether), separated by a lateral tunneling barrier. 
We further consider the situation where there is no cross 
tunneling between the top layer of the system A and the 
bottom layer of the system B. Tunneling occurs either (i) 
between the top and the bottom layer of the same system 
or (ii) between the same layers of the system A and B. A 
schematic diagram is shown in FiglSl 

The configuration similar to the one depicted in Fig|S| 
has been, in fact, proposed previously 3 ^. The important 
twist in this paper is the application of the interlayer 
current flowing vertically through the top and the bot- 
tom layer of one of the bilayer systems. The purpose is 
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Ij = J sin§ 



l=2et<S> sin§ 



FIG. 8: Schematic diagram for the proposed Josephson effect 
in counterfiowing current measurements in a coupled set of 
two bilayer quantum Hall systems. 



to induce a non-zero interlayer phase difference in one 
of the bilayer systems while the other system has none. 
Note that, in bilayer quantum Hall systems, the phase 
of the condensate is controlled via vertical interlayer cur- 
rent flow which affects the internal degree of freedom of 
the condensate. It is interesting that there is no similar 
way of controlling the phase in usual superconductors. 

Specifically, when the applied interlayer current is /, a 
non-zero interlayer phase is induced so that 



sin0 = I/2et(S), 



(26) 



where (S) is defined in Eq. l|18|l . Our prediction then is 
that, due to the analogy with the Josephson effect in 
superconductivity, there will be two counterfiowing cur- 
rents with one flowing between the top layers of the sys- 
tem A and B, and with the other flowing between the 
bottom layers; the Josepson effect for the neutral, total 
current in our system is given by 



I, j = Jsin( 



J 



2et{S) 



(27) 



where the proportionality constant, J, is dependent on 
tunneling parameters between the system A and B. Note 
that, when Ij is the current flowing between the top lay- 
ers, the current of —Ij flows between the bottom layers. 
The net current is zero, but it may be possible to measure 
these two counterfiowing currents individually. 

Let us mention what conditions should be satisfied for 
our prediction to be realized. In the above, we have al- 
ready mentioned that the cross tunneling between differ- 
ent layers of the system A and B is assumed to be neg- 
ligible. This is due to the fact that the cross tunneling 
cannot carry the supercurrent. A maybe more important 
condition, however, is that the applied current, /, should 
not exceed the critical current since, otherwise, the co- 
herence will be lost. The experiment described in the 
above, therefore, can be performed only in a very narrow 
range of the applied current, which, in turn, corresponds 
to a very narrow bias-voltage range. 



VI. CONCLUSION 

It has been shown in the preceding sections that the 
interlayer tunneling conductance peak near zero bias is 



strongly enhanced, but fundamentally finite even at zero 
temperature. The reason why the interlayer conductance 
peak is not intrinsically infinite can be attributed to the 
fact that the experimental setup measuring the inter- 
layer tunneling conductance is not a setup for the true 
Josephson effect in the bilayer exciton condensate. To 
substantiate this idea, we have computed the height of 
the zero-bias interlayer conductance peak as a function 
of interlayer distance, which is, then, compared with ex- 
periment. Considering that the only fitting parameter is 
an overall scale factor which strongly depends on sample 
specifics, the agreement with experiment is reasonable. 
Specifically, our theory provides a reasonably accurate es- 
timate for the critical distance in which the conductance 
peak disappears. Our theory is contrasted from previous 
theories in that it determines the critical distance by ex- 
plicitly computing the interlayer tunneling conductance 
as a function of interlayer distance. Thus, our theory is 
a transport theory for the tunneling conductance, which 
is precisely what is measured experimentally. 

While one of the main goals of this paper is to develop 
a concrete theory of the coherent tunneling conductance, 
it is another goal of this paper to make a sharp distinction 
between the spontaneous coherence due to the Coulomb 
interaction and the induced coherence due to the single- 
particle tunneling. In this respect, the detailed compu- 
tation of the conductance peak height as a function of 
d/ls is also very important because it offers such sharp 
distinction; if the coherence is due to the Coulomb inter- 
action, the conductance peak will completely disappear 
at a finite d/ls, as discussed in SeclTTl 

We have also discussed the similarities and the differ- 
ences between the bilayer quantum Hall interlayer co- 
herence phenomena and the related collective coherent 
physics manifested in magnetic phenomena and in su- 
perconducting Josephson effect. In particular, we have 
proposed an experimental setup for the Josephson effect 
in couterflowing current measurements in a coupled set 
of two bilayer quantum Hall systems, which is a precise 
analogy with the real Josephson effect in superconduc- 
tivity. 
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